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THE CRITICAL OR WHIRLING SPEED. 
OF SHAFTS. 


By Wittram A. NOBLE. 


INTRODUCTION 


TuE problem of the whirling of shafts has received the attention 
of many eminent engineers. Consequently the information avail- 
able is voluminous and the curious find themselves confronted with 
the task of perusing many works. Many of these writings are not 
easily accessible and also contain much information which is of 
purely academic interest. In view of these considerations the 
author has been prompted to write this pamphlet and it is hoped 
that sufficient information has been provided to enable the practical 
draughtsman to approach any critical speed problem with con- 
fidence. 


In the preparation of a pamphlet the tendency is to view the 
subject from the point of view of the industry in which one is 
engaged. In this case, although the author’s experience is with 
the marine steam turbine, it is emphasised that the theory is quite 
general and can be applied to any shaft. Thus those in the aircraft 
and electrical industries will find that the text is applicable to the 
problems encountered in these branches of engineering. 


It is a primary consideration in the design of any power unit, 
be it rotary or reciprocating, that the periodic disturbing forces 
set up when the machine is in operation do not have a frequency 
in the neighbourhood of the natural frequency of the moving parts, 
or alternatively, if the frequencies do approach each other that 
they be out of phase and hence damp out. 


Let it be supposed that a body is vibrating naturally in the 
manner. shown in Fig. 1 (a)—ordinates represent amplitude and 
abscissae time—and in addition to the vibration (a) it is also 
subjected to the forced vibration shown at (b). The resultant 
vibration obtained by adding the ordinates of (a) and (b) alge- 
braically is shown at (c). Since (a) is the natural frequency of 
vibration of the body and (b) is exactly in phase with it, we have 
the condition for resonance and consequently the amplitude of the 
combined vibration (c) will in theory become infinite ; in practice, 
the value would be finite, due principally to the effect of elastic 
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hysteresis in the material of the vibrating body, but the amplitude 


of the vibration may be such as to induce stresses which would 
cause failure. 


Referring to Fig. 2, let (d) represent the natural frequency of 
vibration of a body and (e) a forced vibration exactly out of phase 
with that of (d). The resultant vibration is shown at (f) and the 
damping effect of the vibration (e) is clearly perceivable. 
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Fig. 2. 
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With regard to turbine rotors it is the in phase vibrations 
(Fig. 1) which must be guarded against, and it will be appreciated 
that it is imperative that the rotor must be run at such a speed 
that the forced vibration set up by the eccentticity of the mass 
centres does not approach the critical or natural frequency of 
vibration of the shaft. In reaction turbines, the small clearances 
of both axial and radial clearance blading indeed precludes the 
existence of any appreciable vibration and in view of these con- 
siderations it is the practice to run the rotors of large turbines 
somewhere in the region of 30% to 50% below the critical value. 


The critical speed of a turbine rotor may be defined from con- 
sideration of the following— : 


In any shaft, be it in a horizontal or vertical position, it is 
impossible to have the centres of mass of the parts situated in the 
axis of rotation. In the case of the horizontal shaft, the bending 
due to the weight of the shaft prevents the attainment of this 
condition and in the case of the vertical shaft there is never exact 
coincidence between the mass centres and the axis of rotation, 
due to inaccuracies which occur in the course of manufacture and 
which are incapable of being eliminated completely. As the shaft 
speeds up from rest, these eccentricities cause a centrifugal force to 
act on the shaft causing it to bend and ultimately a speed is attained 
at which the centrifugal force is equal to the elastic righting force 
exerted by the shaft. If this speed is maintained, the shaft will 
fail and it is said to have failed at the fundamental critical speed. 

‘ 

The unstable condition just described requires time to develop 
and if the shaft is rapidly accelerated through the fundamental 
critical speed petiod it will steady up again and can be run with 
safety. In small turbines, that of de Laval for example, this is 


the practice and in Section 9, higher critical speed values will be : 


discussed briefly. 


SIMPLE THEORY 


The simplest case is that of a light shaft with a comparatively 
heavy wheel, as depicted in Fig. 3. 


Assuming the shaft to’ be weightless, 


Let W = weight of wheel—lbs. 
= angular velocity—radians/sec. 
deflection at wheel—inches. 
lateral stiffness of shaft—Ibs./inch deflection. 
= gravitational acceleration = 386-4 inches/sec.?. 
‘ B 


Rr & 
iT 


“ 
t 
3, 

3 


‘ 


The eccentricity is y and the centrifugal force acting on the 
shaft is given by : 


ge 
, ,a oe a a) 7 
The elastic righting force exerted by the shaft is given by : 
Fz, = P.y (2) + | 
The condition for the critical speed is that F. = Fy. ty * 
Woaty=P.y iy 
& 


we = critical speed—radians/sec. 
Ne = critical speed—revs./min. 


‘ 1 ‘ : i 2 a 
Ea ee eee eee eee een Pole ee 
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In order to evaluate ic, it is necessary to determine‘the value 
of P. 


Referring to Fig. 3, we proceed to evaluate P as follows : 


The bending moment on the portion of the shaft to the left of 
W is, for any distance x from C as origin, 


M. = B.W.* 
Mr = — L 
and for the portion to the right of W is, 
B.W. 
M.= - Lot We -A) 


The equation for the deflection is, 
r M 

y | —— dx . de 

. JJ. I 


= D* for solid shafts and 
64 


i 


where I 


a on 


a (D*~d*) for hollow shafts, 


j D and d being the shaft outside and inside diameters respectively. 


’ For the left-hand portion of the shaft, the deflection equation 


expands to 
B.W.x# 
we Fae 
and for the right-hand portion is 
BLW. x8 2 W x8 W.Ax? H 
ae” Bite eh. OL. PP 


The four constants I, G, H and K are found from a knowledge 
of the boundary conditions. 


Assuming the shaft to be simply supported, the boundary , 
conditions are : 

(a) The deflection is zero at C, i.e. when x = O. 

(b) The deflection is zero at D, i.e. when x = L 


(c) The deflections for the two portions are equal at N, 
te. when x = A, 


(d) The slope of the elastic curve is the same at N for the 
two portions, 7.e. when x = A. 


a ‘ 


aa 
% 
Wa 


. 
a” 
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The values of the constants are : 
Ww. aa WoL L wat" 
F= 7 (- S45 A) “21.E 
6L.I.E 2.1.E 3 3 is 


G =0  % , 
¥, 3 J 
Hy . WA _WLL (-2+4-a) 
6L.1.E 21.E 3 3 
E Ww. As : ‘A 
K = = > i 
6.1.E 
Hence the deflection’ for the left-hand portion is given by 
oo | Bx (L?- x? BY 6) 
61.E 1 
and for the right-hand portion is : 
Ww ( B.x (L?- +*- B?) ) 
=o COS —A)3 6 
de, 61.E i L + (x ) j (8) P 
The deflection at N by either (5) or (6) is : ; 
W. A*®. B® 
= 7 
Ye = 3.1.E.L> vi ‘ 
When yx’ = 1, W=P ; 
32T.s Bok 
= Be 8 
e A®. B? 8) P : 
Hence from (4), == ‘ 
2...E.1. : 
iene =k an) 


7 © W.A2B 
It may be easily shown that the critical speed is numerically 
equal to the natural frequency of vibration of the loaded shaft. 


Referring again to Fig. 3, let it be supposed that the shaft is 
placed in a vertical position and let N be displaced a distance y 
horizontally and then released. The shaft will then vibrate with 
simple harmonic motion, and the amplitude of the (assumed un- ‘9 
damped) vibration will be y at N. 

Let 2 = frequency—oscillations/sec. 


The kinetic energy of the wheel in the mean position is / 


Ww 
KE. = — (2 7.y.n)° 


At the extreme position, the kinetic energy is converted into 
strain energy. i. 


The strain energy is } P. y? 


. 
. 


oo ee Pe 
tN el 
rot 
i 
s 
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Equating the strain energy to the kinetic energy, 
Ww 
4 P.y?= — (2 7.9.0)? 


2g 7 
1 ya 
We 5 WwW 


T 


Let nu. = frequency—oscillations/min. 

2 Pig 

Then — 30 — 
nF W 

Substituting the value of P from (8), 


_ 30 «/81-E.L.g 


(10) 


Comparison of (9) and (10) proves the identity. 


Example—Estimate the critical speed of a shaft $” dia. and 27” 
long when it carries a load of 30 lbs. at a distance of 9” from one. 
end. Neglect the weight of the shaft and assume it is simply 
supported on end bearings. : 


The values of the various quantities are : 


I = -0155 inch! E = 30 x 108 Ibs./sq. inch. 
L = 27° g = 386-4 inches/sec.? 
W = 30 lbs. A = 18” 

B= 9° 


By formula (9) the critical speed is, 
30 «f 3 x 0155 x 30 x 10° x 386-4 x 27 
a ee 30x324xSl 

No = 1,295 RPM. 


It should be noted that although Ib./inch units have been used 
in the foregoing, the theory is quite general and equally applicable 
to the metric system of units. Lb./inch units are used throughout 
the work, but the reader will find no difficulty in inserting the 
appropriate metric units should the need arise. 


DUNKERLEY’S EMPIRICAL FORMULA.* 


In the case of a shaft and disc rotor with a uniform shaft, an 
empirical method due to Professor Dunkerley, may be employed 
to calculate the critical speed. This method has been confirmed 
by experiment and may be applied safely to a two-bearing rotor 
supported on end bearings. | It should, however, never be applied 
to a multi-bearing rotor as the results are liable to be very inaccurate. 


* For a theoretical explanation of Dunkerley’s formula, see “‘An Analytical 
Method,” page 18. 
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If the shaft is not uniform, a method to be described later should 
be employed (page 12). : 
Let N = critical speed of whole rotor—revs./min. 
n, = critical speed of shaft alone—revs./min. 
Ng, Nz, etc. = critical speed due to each wheel mounted sup- 
: posedly on a weightless shaft having the same 
lateral stiffness (the same dimensions) as the 
actual shaft—revs./min. 
1 1 1 1 
Ne ny y Ne tg? Fea G2) 
The critical speed due to each wheel is found by the method 
already described. The critical speed of the shaft alone is found 
as follows : , 


Then 


- Let w = weight of shaft—lbs./inch length. 
E = Young’s modulus—30 x 108 Ibs./sq. inch. 
w = angular velocity at critical speed—rads./sec. 
y = deflection at any distance x from a fixed origin—ins. 
1 = length of shaft between supports—inches. 
§ = gravitational acceleration = 386-4 inches/sec.? 
I. = moment of inertia of shaft cross-section—inch.4 


ny critical speed—revs./min. 
The centrifugal force tending to bend the shaft is, 
w 


Fy = = w*. y per unit length. 
The bending force is also given by 
dy F 
F, = E.I —— per unit length. 
dx! 
4 2 , 
‘ Hence fA 
dx4 g. Et 
Let mo = 
: RirakSve DT 
, ae 
Then Ss why (12) 
dx4 = 
enares 
and o = A, ob (13) 
w 
In order to evaluate m2, it is necessary to partially solve equa- 
tion (12). 


Substituting y = A . e®*, where ¢ zs hyperbolic logarithm base, 
the auxiliary equation is, 
#14 = a 
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The roots of this equation are, @ = + m°and @ = +1%m, 
where i = V -1. : 
oy = Ae™ + Bo e-™™ 4 C. emi+ Demmi 
| This may be written, 
y =F sinh mx + G cosh mx + Hsin mx + J cos mx (14) 


The boundary conditions for a simply supported shaft are, 
(a) The deflection is zero at the supports, 1.e. when x=0 and 


. et. 
; (b) The bending moment is zero at the supports and hence 
oe, = 0 when x=0 and x=1. ! 
dx® 
i When x=0 in (14), G+J=0 
Differentiating (14) twice and substituting x=0, G - J =0 


Hence, G=J=0. 
When x=/ in (14), F sinh ml + H sin ml = 0 
Differentiating (14) twice and substituting « = L 
F sinh ml - H sin ml = 0 : 

Hence, 2H sin ml = 0 

If H is not zero, the shaft is vibrating and sin ml = 0. 

ml = 7, 27, 3 7, etc. 

The value mi = 0 is not admissable as the shaft would then 

be at rest. 


Hence, mm = ae ae om etc. 
Jt L L 
Substituting the values of m in (13), it will be seen that there 
are an infinite number of critical values and they are proportional , 
to the numbers 1, 22, 3°, 4%, 5?, etc. 
“Selecting the lowest value, the fundamental critical speed is, 


} . 7 g.E.D 
b eee V w Gs) 


30 7 VJ g.E.1 7 
: and nm = a ae © (16) 
p) Example—Estimate the critical speed of a shaft 2” dia, and_ 


48” long when it carries a load of 40 lbs. at a distance of 15” from 
one end. Assume that the shaft is simply supported on end 
bearings and take the weight of a 2” dia. shaft to be 0-9 lbs. per- 
inch length. 

The values of the various quantities are, 


I = -7854 inch‘. E = 30x 10° Ibs./sq. inch. 
L = 48’ g = 386-4 inches/sec.* 

W = 40 lbs. 5 A = 15’ 

B = 33” w. = -9 Ibs./inch. 


-. 
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By formula (f6) the square of the critical speed of the shaft 
alone is, - 
, _ 900 a x 386-4 x 30 x 10°x-7854_ 
te 48 x 9 


16962138-7 


By formula (9) the square of the critical speed of the load alone 
supposed mounted on a weightless shaft is, 
900 x 3 x +7854 x 30 x 108 x 48 x 386-4 
Me a2 x 40 x 295 x 1089 


i 


= 12197384-9 


By formula (11) the critical speed N of the whole shaft is found 
from 


1 1 1 


= + = 
Ne 16962138-7 12197384-9 
and N = 2660 R.P.M. 


THE TWO-BEARING ROTOR. 


It has been pointed out by Lord Rayleigh that within wide 
limits different assumptions as to the shape of the deflection curve 
for a rod making free transverse vibrations do not materially 
affect the calculated value of the natural-frequency. The correct 
value is a minimum and all others are in excess of this. 


Advantage is taken of this fact in calculating the critical speed 
of turbine rotors. If the curve of static deflection is assumed to 
be the shape of the actual deflection curve when the shaft is at 
speed and the total resilience of the shaft is equated to the total 
kinetic energy in the mean position a close approximation to the 
actual critical speed will result. In the case of the two-bearing 
rotor, simply supported at its ends, the approximation is within 
1%, of the actual value and for all practical purposes this is suffi- 
ciently accurate. 


Let w = variable load at any point along the shaft—lbs. 
y = deflection at any point along the shaft—inches. 
wc = critical speed—radians/sec. 

Ne = critical speed—revs./min. 
g = gravitational acceleration = 386-4 inches/sec.? 

The resilience of the shaft at any point is, 

y=hwiy 

The total resilience is, 

R=42w0.y (17) 


, * 
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The kinetic energy for the mean position at any point is, 


k= —.are. 3? 
2g 
The total kinetic energy is, 
K= 2 sw.y? . (18) 
2¢ . 
Equating (17) and (18), 
we? 
1l2w0.y = aw. 
2g 
[ZSw.y.g 
t f= VS 19 
then, a, Vv Sw. (19) 
30 Zw.y.g 
and jig = = —— = , (20) 


7 aw. 


The problem of a complex shaft may be treated graphically. 
A simple example of the method to be employed to find the curve 
of static deflection will suffice to render the reasoning apparent. 


In Fig. 4, ab represents the line of action of a force F. It is 
required to determine the moment of F about a point c. 


L Fi 
= 4 
y 


Fig. 4. Cc 


a 
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Draw de parallel to a b and along d e, measure a distance g h to 
represent the magnitude of F. Choose any pole & and join kg and 
kh. Through c draw /m parallel to ab. From any point x in ab 
draw 2p and xq parallel to kg and kh respectively. Then gf re- 
presents the moment of F about ¢ if a suitable scale is employed 
for measuring its magnitude. 


Proor.—Since the triangles 7 pq and & gh are similar, 


t s 
Pg gh 
where ¢ and s are the altitudes of the two triangles. 
Hence, t.gh = s.pq (21) 


Since gh represents the magnitude of F, then ¢ . gh is the moment 
of F about c, and if s=1, 
= pq 


Thus with a suitable scale, pq is the moment of F about c. 


ScALE.—Let s=x times the linear unit employed. Let ¢ be 
drawn to a scale where one linear unit represents y linear units 
and gh be drawn to a scale where one linear unit represents z units 
of force. 


By (21) t.gh=s.pq 
Applying the above ane 
t.y.z.gh | = - PY 


Hence the scale of pg is ¥ ¥ a4 Panis of moment per linear unit. 


The foregoing construction is quite general and may be applied 
to the case of a complex shaft. The inclusion of such elementary 
reasoning will be justified if it eliminates the errors which often 
occur in arriving at the scale of the bending moment diagram. 

The change of slope in radians between any two points on a 


beam of constant section is given by, 
XX 


ay = = | M. dx 
dx Bill dy 
where M = bending moment. 
E = Young’s modulus. 
I = moment of inertia of the section. 


.% = distance from a fixed origin. 
= deflection at distance +. 


Xy 


The quantity | M . dx represents the area of the bending 


moment diagram between the points considered. 
Actual area of B.M. diagram x constant 
Ek 
gives the change of slope in radians between the points considered. 


Hence, 
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The constant is found from a knowledge of the scales employed 
in drawing the bending moment diagram. 


Let 1” vy inches ~ length. 
iL z lbs. — force. 
Let the polar distance = x inches. 


The scale of the ordinates of the bending moment diagram is 
as has been shown previously : 

1” = x.y. 2 1b. inch units. 

Hence one square inch of the bending moment diagram re- 
presents, 

y (%.a.2) =x y? 2 Ib. inch® units = constant. 

If the slope values are considered as imaginary loads acting 
through the centroids of the bending moment areas, the deflection 
diagram may be drawn in a similar manner to the bending moment 
diagram. ; 

Referring to Fig. 4, let gh represent the change of slope between 
any two adjacent points on a beam. Then gp measured with a 
suitable scale represents the change of deflection between the 
points considered. 

Let 1” y inches — length. 
1” z radians — slope. 
Let the polar distance = x inches. 
Then the scale of the deflection diagram is 
1" = x.y. 2 inches. 


tot 


@ a & $ 


Fig. 5, 


The bending moment areas with the exception of the triangular 
areas at the extreme ends approximate to a trapezoid. Fig. 5 
illustrates a useful graphical method of determining the centroid 
of a trapezoid. 


Let abcd represent any trapezoid. Produce ab and ed to 
e,f,gandh. Make ae and bf equal to ed and dh and cg equal to ab. 
Join eg and fh; kis the centroid of a bed. 

In order to illustrate the method, an example has been worked 


in full and the procedure may be summarised in the following 
manner. Reference is made to Fig. 6 and Table 1. 


* eS mm be Adee 


ee ee 
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(a) Draw the shaft to as large a scale as practicable. (In the 
example the scale is 1"= 1 ft., the scale being kept small to facilitate 
reproduction). 

(b) Divide the shaft into a number of sections, I having a 
constant value over each section and estimate the weight of each ; 
blading, etc., must be included. Let these weights act through 
the centres of mass of the sections. 

(c) Draw the force polygon, using as large a force scale as 
practicable. In Fig. 6, the scale is 1”=1000 lbs. and the polar 
distance is 3”. "i 

(d) Project the load lines vertically downwards, and by drawing 
lines parallel to the rays of the force polygon to intercept the load 
verticals, construct the bending moment diagram. 

(e) Estimate the areas of the sections of the bending moment 
diagram under each load (sq. ins.). 

(f) Calculate the change of slope over each section from the 
equation, 

i Actual area of B.M. diagram x constant 
G2 . 

The value of the constant is, 

3 x 12?x 1000 = 432,000 inch? lbs. 

(g) Using the slope values, draw the imaginary force polygon. 
The scale here is, 1”= 1/50000 radian. 

(h) Draw verticals through the centroids of the bending 
moment areas and from the imaginary force polygon construct the 
deflection diagram. 

The scale of the deflection diagram is : 

1”= 3 x 12 x 1/50000 = -00072 inch. 

(j) Estimate the deflection at each load and calculate the 
values of w. yandw.y?. The various values are given in Table 1. 

(k) Calculate the critical speed from formula (20). In the 
example, this is : 


30 me g 
tm =— ————— 
7 -0031 


= 6750 R.P.M. 


It should be noted that this is a first approximation, but, as 
has been mentioned already, is sufficiently accurate for a shaft 5 
of this type. 

A more accurate value may be obtained by drawing a new 
deflection curve, using the original deflection curve as an assumed. 
one, calculating the centrifugal force at each load about its respec- 
tive radius of rotation (the respective deflection y on the assumed 


ote, JARRE eee 
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deflection curve) and using these centrifugal forces instead of the 
origina] static loads. Formula (20) then takes the form, 


30 
nm = — 


(20a) 


where P 


centrifugal loading at the speed ic. 


= deflection on the derived deflection curve due to 
the load P. 


AN ANALYTICAL METHOD.* 


The semi-graphical method of calculation just described will 
undoubtedly appeal to many, but some may prefer an analytical 
method. In any case an alternative method provides a ready 
means of checking the result obtained by the graphical method 
and also the method about to be described is not so laborious. 


The problem of determining the critical frequency reduces to 
that of finding the deflection at various points of a beam when it is 
subjected to some particular system of loading. The system of 
loading as far as the critical speed calculation is concerned is 
dependent on the manner in which the shaft is supported. 


WW. Ws Wy 


L 


Fig. 7. 


Consider Fig. 7. 
Let yx = deflection at load Ws. 


Then the work performed by Ws when the beam is deflected a 
distance yx at this point is Ws/2 . yy. 


* “The Critical Speed of a Turbine Rotor,” by J. Morris and S. J. E. Moyes: 
Engineeving, 30th June, 1939. 
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Let Wx be increased to (Wy +1). The work done will then 
be § (Ws +1) (yx +”), Where y is the deflection due to the unit load 
acting alone in the position of Ws. 


Ws . Vs Ws. ¥ ayy y 
+ + 


i 


3 (Ww +1) (ye +9) 


2 2 2 2 
but y = s/Ws by 
ys = Wye y * 
Hence the work done is q 
Wr. yx y ‘ 
5) + yy + i (21) | 
This work is stored as energy by the shaft. 
Let M = bending moment at any point of the shaft due to 
the normal loading Ws. 
m = bending moment at any point of the shaft due to 
unit load acting alone at the position of Ws. 
E = Young’s modulus. 
I = variable moment of inertia of the shaft cross- 
section. ? 


The change of slope of the deflection curve for the element @x is, 
(M_+ m) 
———_ 
Eon 
The energy stored due to this change of slope is, 
(M + mt)? 
2Ew 


\ 
The total energy stored due to the load (Ww +, 1) acting through 
the distance (ys + 9) is, 


dx 


* (M+ mm)? 
P ?e.0 | 
2 [ M ee [ M.m ae rn rene iy 
~ | 2E.1 © i eT J,2E.1 i 
Consideration will show that, 
Wy» Vs e Me 
, £ [ 2E,1 
2x — [ae 
“iii g °° | oe ™ (22) 
u M.m 
hence ys = EI dx (23) 


= > 


ia _ =. 
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If the sign of integration is replaced by the symbol & in (22) 
-or (23) we have the approximations : 


Som 


i oe ax By 

’ >. Bal a) 
s M.m 

te = BS ae (25) 


where N is the number of positions considered. 


' Thus yx can be calculated to any required degree of accuracy 
by splitting the shaft into short lengths and calculating the mean 
values of M, I and m for each length. The deflections thus found 
may be substituted in formula (20) and the critical speed calculated 
as before. 


In order to simplify the work, a formula for the critical speed 
was derived from the rational explanation of Dunkerley’s empirical 
formula. 


The method of explanation of Dunkerley’s formula used here’ 
is due to Hahn* and in order to render the reasoning apparent a 
simple example is taken ; the result, however, is quite general. 


Ww, We 


Fig. 8. 


Let any shaft AB, Fig. 8, be loaded, as shown, and let the 
critical angular velocity of the shaft be wx. 


Let the deflection at the critical speed in way of each load 
be y, and yp. The centrifugal loading then is, 
‘4 
ok - we? . yy, and 
g & 
Let the shaft be at rest and let all the loading be removed. 


. Wy + Me 


‘Then unit load in the position W, will produce the deflections 8y1 


and 8,, and unit load in the position W, will produce the deflections 
So, and Soo. 


* «Steam and Gas Turbines.” Stodola. Page 1110. 
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Hence the deflections at the critical speed are given by, 
W, Fi We 2 
dy = — on? dy fe =. an? do. - Vy 
& g 
W. A Vo 2 
vay = na Bia Nyt =. wx". 329. V2 


These may be written 


W. 2 2 
( : an? By ~1) 94 a Vs an? . Bo. Ve = 0 
g g 


Ww, W, 


ar? . Bo. Vy + ( 
zg 
56 


, The eliminant of the equations for finite values of y, and Vo 
is given by the determinant, 


g 
é 


+ WK? . B21) da 


WwW. y 
( + ay? By.-1 ) + We cid By, 
§ g 
= 0. 
y » 
Me oat dag (ME seat tee 1°) 
§ : g 
Expanding we have, 
W W. W. : 
i’ » wx! 84. 8y2- z wo? . by + 1 
We Vy 
- —2 2. Sao Wes : MY ile Ba Rigen B 
g g g 
Approximating to the first order of small quantities, 
1 W, WwW, 
ee ‘ —  . Sy. > 2 
, aa 8 + z Soe (26) 


If in the example we consider the critical angular velocity for 
the loads acting singly, then for load Wj, 


Ws 
(Sheet bn - 1) va = 0 


W, 5 
7" w,.5, -1= 0 

hence, E mes: 8 (27) 
wy & 
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d Similarly, : 
e _ We | Bon (28) 
ws g§ 
Substituting the values of (27) and (28) in (26) we have, 
1 1. 1 
3 oy Peg (29) 
wy” a," We 


which is the equation of Dunkerley. 
Equation (26) may be written, 


ok a Bee 
4 wie g 
, i, © a 
‘ “: 2W.8 
Let ue = Critical speed — revs./min. 
5 oS 
7 ZW.8 iF) 
The values of 8 can be calculated from 
\ Sn 
S=y = 2 eau (24) 
Equation (24) can be written in a form which lends itself to 
‘ systematic tabulation. The bending moment for unit load at Z, 
Gz. ba 


' Fig. 9 for a beam of length L is at Z. 


The bending moment for all points to the left of Z is given by 
a.,/L and for all points to the right by az. b/L. Inserting these 
values in (24) and taking due note of the limits we have, 

b2 © ia? det a? Nb. dx 
& = z E 
2 Ee “1 1) * Ee ad 


Fig. 9. 
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Further, 
Wa. be 2 a®. dx Wade? 8 2 dx 
85 SS : 9 
We Bae ee E.L? ‘sa: I (1) 
Also formula (25) may be written, , , 
bz se M.a.dx > s M.0b.dx 


The critical speed of the shaft of Fig. 10 has been calculated 
by formula (30) using the values of (31) and also by formula (20) 
using the values of (32). The values of the various quantities 
for the two methods are given in Tables 2 and 3 respectively. 

The critical speed by (30) is 17,170 R.p.m. and by (20) is 17,680 
R.p.M. Assuming the latter value. to be exact, that given by (30) 
is about 3% too low. The value given by (30) is always on the 
low side and may be as much as 5% removed from the actual value. 


THE OVERHUNG SHAFT. 


The overhung shaft is seldom encountered in marine turbines, 
but it is met with in aircraft gas turbines and certain types of 
electrical generating plant, such as water wheel driven alternators. 

The procedure is essentially the same as that for the two-bearing 
rotor, but there are a few points which require special mention. 

Referring to Fig. 11, it will be seen that the deflection curve 
dips between the supports and rises outside the right-hand support 
at a particular instant when the shaft is at speed. It is important 
to note that this is the condition which gives the lowest or funda- 
imental critical speed, that is, a two-node deflection curve for the 
shaft, and it will be appreciated that in order to attain’ this the 
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loading inside and outside the supports must be assumed to act in 
opposite directions. Hence if the loads inside the supports be 
considered positive acting downwards, then those loads on the 
overhung portion must be considered negative acting upwards. 
It will thus be seen that it is necessary to have some notion of the 
shape of the deflection curve when the shaft is at speed and this 
is the case for all rotors except the simple two-bearing one of 
Fig. 6, where the loading is all of the same sense. 


The force polygon is drawn as before, due note being taken of 
the signs of the loads. The bending moment diagram is constructed 
by taking the rays of the force polygon in order and drawing 
parallels to intercept the load lines. In the first instance the 
right-hand support line must be ignored and when the extreme 
right-hand load line is reached it will be found that the ray 9 is 
left over ; a line must be drawn parallel to this from the extreme 
right-hand load line to intercept the right-hand support line. The 
diagram is closed by joining AB. 


The deflection diagram is drawn as befcre, the only difference 
being that the diagram is closed by the line CD drawn through the 
points of zero deflection. 


The critical speed of the shaft of Fig. 11 has been worked out 
and the values of the various quantities are given in Table 4. The 
critical speed is by formula (20) 11,085 x.p.m. If any doubt should 
exist regarding the accuracy of the result the work will have to be 
repeated, using formula (204). In general, when the overhung 
load is relatively heavy the latter calculation will have to be 
performed. : 


THE MULTI-BEARING ROTOR. 


It is necessary before embarking on a discussion of the method 
to be employed in this case to state exactly what constitutes a 
multi-bearing rotor and what does not. A multi-bearing rotor 
may be defined as a continuous shaft or a series cf shafts rigidly 
coupled together and supported on three or more bearings. If 
flexible couplings interpose the lengths between them must be 


considered as separate units for the purpose ef the critical speed 
calculation. . 


The problem of the multi-bearing rotor reduces to that of 
finding the deflections at various points of a propped beam. The 
method of determining the deflections used here is that due to 
Wilson.* 

Consider the beam ABCD, Fig. 12. Let the supports B and C 
be removed. The beam will deflect, as shown at (a) and the 
deflection at B is ys and at C is ye. Now, let all the loading be 

* Proc. Roy. Soc., vol. 62, Nov., 1897. 
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Fig. 12, 


removed and unit load alone applied at B will produce the de- 


flections 8,, and 6,, at B and C respectively (h). 


process by applying the unit 
at B and 64, at C (c). 


Let the reaction at the sup’ 


since the deflection of a beam 
Vo 

Me R,. 

The above equations can | 
are # intermediate supports t 


Let the beam be loaded at 


The deflect 
deflections at a are considere 


R, respectively. 


and they are added algebrai 
of the beam will result e. I 
diagram d can be constructed 
and the reactions R, and Ry. 


In Fig. 12 the loads hay 


R,. 8; 
S12 


Repeat the 
oad at C and the deflections are 84, 


port B be R, and at C be Ry. 
is proportional to the load, 
+ R..8y) 

+ Re. doo 

be solved for Ry and Rg, and if there 
here will be 1 such equations. 

B and C with loads equal to R, and 
ion diagram is shown at d. If the 
d positive and those at d negative 
cally, the actual deflection diagram 
should be noted that the deflection 
from the deflection diagrams b and c 


Then, 


1 


e all the same sense, whilst for the 


purpose of calculating the critical speed the loading is considered 


~<a 
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Arrows Dewors Osrectian Or Loaoine 
Fig. 13. 


alternately positive and negative over the consecutive spans from 
consideration of the shape of the deflection curve at speed. 


The shape of the deflection curve and the sense of the. loading 
for a three, four and five bearing rotor is shown in Fig. 18. 

It is unnecessary here to work an example, but the method may 
be summarised as follows :— 


(a) 


(g) 
(h) 


Remove all the supports except the end ones and give the 
loads signs in conformity with Fig. 13. 
Draw the deflection diagram in the manner which has 
been described for the two-bearing rotor. 
Draw deflection diagrams for unit load acting alone at 
each of the removed support positions. 
From the general equation, 

y = R,.8y + Re. dg +... + Re. yy 
form N equations and solve for Ry, Rg... Ry. 
Construct the deflection diagram for the loads R,, Rg, 
. .. Ry acting on the shaft. 
Draw the actual deflection diagram for the supported 
shaft by subtracting that of e from that of 8. 
Calculate the values of W . vy and W . 4°. 
By formula (20), calculate the critical speed of the shaft. 


The value calculated by (20) must be considered in this case 
as a first approximation, and it is therefore necessary to repeat 
the work using formula (20a). If the derived critical speed is not 


“= 
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reasonably near to the original value (say 5% difference) the 
calculation will have to be performed a third time. It should be 
noted, however, that in practice it is seldom necessary to perform 
this third calculation. : 


THE GYROSCOPIC EFFECT OF THE ROTATING MASSES. 


In the discussion up to the present it has been assumed that the 
masses have been concentrated at points situated in the axis of 
symmetry of the shaft ; it is apparent that this is not the case. 

When the shaft is at speed a mass, such as a wheel, is inclined 
at an angle to the axis of rotation, the inclination depending on 
the slope of the deflection curve and consequently a moment is 
exerted by the wheel due to the tendency which it has to rotate 
in a plane perpendicular to the axis of rotation. This moment, 
according to its sign, tends to increase or decrease the deflection 
of the shaft and hence lower or raise the critical speed. 

Cases may occur where the moment is appreciable; a heavy 
wheel close to the bearing for example, or a wheel at the end of an 
overhung shaft, and it is therefore necessary to show how its effect 
may be taken into account. 

In general, the gyroscopic moment has a decreasing effect on 
the deflection of the shaft and it is proposed herein to confine the 
discussion to this aspect of the problem. . 


ole 


dm odx six 2 


. Save F 
IMOX CoB 


Fig. 14. 


The wheel, shaft and deflection curve all rotate with the same 
sense and with the same angular velocity w (Fig. 14), and conse- 
quently only the centrifugal forces are to be considered in calcu- 
lating the gyroscopic moment. 

. The inclination of the wheel to the vertical is a radians and 
since this is a small quantity we may write, 
cosa = 1 andsing = 0. 
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The centrifugal force for the element of mass dm is 6. wx. * \ 
The vertical component of the centrifugal force is, 
SF y = dm. w%.« cos B = Sm. w* (y +2) 
and the horizontal component is, 
8 Fy = du. w% x sin B = dm. ws 
The total vertical force is, 
Fy= 28m. w? (yt2) = Zim. wry + Zbm. wz 
but, 2 dm .w*. z = 0, since the simple moment of mass referred 
to a diameter passing through the centre of gravity of the wheel 
is zero. 
The total vertical force is, therefore, 
Fy = m.wiy 
The total horizontal force is, 
Fy, = 2 6m. w%s = 0. 
P The moment of the centrifugal forces about a diameter situated 
in the centre of gravity plane of the wheel is, 
M S8Fy.z.a 
ZS dm.wiy.za + ZF bm.0% 2a 
Z bm. w% za 
The quantity 2 5m. 2? is approximately the mass moment 


of inertia of the whcel referred to a diameter situated in the centre 
of.gravity plane. 


ages 
‘ Wee Wie 


i 
cu 


a 


A 
rot 


Hence, M =wia.Ip (33) 
where, Ip = mass moment of inertia of the wheel referred to a . 
diameter passing through the centre of gravity of ‘ 

i the wheel. 
The foregoing reasoning holds good for any of the rotating 
masses, be they wheels or portions of the drum. 4 


It is a simple, if somewhat tedious, matter to take this moment 
into account. Let it be supposed that in the first instance the 
critical speed has been calculated, neglecting the gyroscopic moment. 
If the value obtained and the derived deflection diagram are taken 
as a basis for calculation, the gyroscopic moments can be computed. 


Fig. 15 shows a two-bearing rotor having three wheels which 
exert moments M,, M, and Mj. The length of the shaft is L. The 
force which must be exerted at the bearings in order to counteract 
the moments is given by 


f, 


, 


Re 2M 


where 2 M is the algebraic sum of the moments M,, M, and M3. 
The ordinates of the bending moment diagram are found in the 
usual manner by taking the algebraic sum of the moments to the 
tight or left of any given section. The sudden changes of the 
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ordinates occur at the points of application of the moments M,, | 
M, and Mg, and it is therefore necessary when calculating the | 
“ ordinates to first of all calculate the moment at each wheel, ne- 
glecting the moment exerted by the wheel and then to repeat the 
calculation taking the wheel moment into account. 

The next step is to correct the bending moment diagram which 
has already been drawn in connection with the calculation for the 
critical speed when neglecting the gyroscopic moments by sub- | 
tracting the diagram for the moments from it. The corrected 
bending moment diagram is then used to derive a new deflection 
diagram and the critical speed is calculated by formula (20). It 
will in general be unnecessary to repeat the calculation, using | 
formula (20) unless the gyroscopic effect is considerable. 

_ The method of correction for a multi-bearing rotor is very 
similar to that just described and is briefly as follows :— 

The critical speed is calculated neglecting the gyroscopic 
moments and then using the results as a basis of calculation, the 
gyroscopic moments are calculated. 

Let it be supposed that the intermediate bearing supports are 
removed and a bending moment diagram for the gyroscopic moments 

. (Fig. 16) constructed in the same manner as for the two-bearing 
rotor (Fig. 15). 
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Fig. 16. 


The bending moment diagram which has been drawn for the 
whole shaft when neglecting the gyroscopic moments and when the 
intermediate bearings are supposedly removed is corrected by 
subtracting that for the gyroscopic moments from it. The actual 
deflection diagram for the shaft is then derived by Wilson’s method. 
The critical speed is calculated by formula (20a) 


CRITICAL SPEEDS OTHER THAN THE FUNDAMENTAL. 


Critical speeds other than the fundamental are seldom required 
except in connection with small high-speed turbines and certain 
categories of electrical generating plant, in which it is desirable 
and practicable to make the rotor as light as possible. In large 
installations, such as are used for ship propulsion, the rotors are 
run at speeds well below the fundamental critical value. 


The actual shape of the elastic curve (the position of the nodes) 
for the higher critical values is a matter for conjecture in all but 
the simplest cases and hence formulae (20) and (20a), which give 
progressive approximations for the fundamental, are no longer 
convergent when applied to the higher categories. 


It is necessary to have some conception of the shape of the 
shaft for the higher frequencies and Fig. 17 illustrates diagramma- 
‘tically the probable shape of the elastic curve for a two-bearing 
rotor when running at the first, second and third critical speeds. 
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Narows Nenote Migection OF Loasing. 
Fig. 17. 


ri “ The method used here for computing the higher frequencies 


~ was suggested by Lord Rayleigh in his ‘Theory of Sound,” and 
5 although the procedure will be found tedious, it gives very reliable 
A results, . 
4 Let w, = first critical angular velocity. 
ws, = second critical angular velocity. 
| Ji = deflection at any point of the shaft at the first critical 
speed. 
; My = Adflection at any point of the shaft at the second 
= critical speed. 
r. w = weight of shaft per unit length. 
. E = Young’s modulus. 
Bi I= moment of inertia of the shaft cross-section (bending). 
4 & = gravitational acceleration. 
79 x = distances from a fixed origin. 
4 L_ = length of shaft. 


Considering the first critical speed, the condition for indifferent 
equilibrium is given by the differential equation, 


d4 
El [3 = = wy? Vy (34) 
56 


dx* 


Similarly, the condition for indifferent equilibrium at the, 
second critical speed is, 


» " fF - ee , “ "i oS Be 
‘ a “ es had ay 


CRITICAL OR WHIRLING SPEED OF SHAFTS BS 
a) y, ew 

i El [22 == w,2. yo 
. ee (35) 

7 Multiply (34) by y. and subtract (35) multiplied by y,; the 

whole is then multiplied by dx and the result summed over the 
length L of the shaft. 
f The result is, 
mo . dt y L . ayy “ 
t ‘ a Tat Reeeltt - few axt ot ci : 
& Ly egy) : (36) 
—.w,?. year | — . We. Vo. Vy aX 
ty L a a V1 Ne .8 a Va Va ) 
2. Considering the left-hand side of (36) and taking the integrals 
Ne . one at a time, we have, integrating by parts, 
x | “eM oy, ad E @y; J 
- E... Ng. dx=]EI. =~. vy = 
( ; ant Yo. 4% ae Me 

2 dey . a} ps fr I yy, Pye dx [a2 

: dx® dx «| a dx” 

2 Similarly the second integral expands to 

L qi Vo By, it 
‘ - Ba. iat 13. de = [ie Te “4 |? 

: @y, dy, T- fea @y, dy, ao (38) 
dx? dx «| aR det . 
mT Inserting the values of (37) and (38) in (36), we have, since the 

quantities inside the brackets are zero from consideration of the 
7 boundary conditions for a freely supported shaft, 
"74 lw Lowy ‘ 
‘ [i Set a sede - [i Fete nae = 0 
i L 
| iti de = © (39) 
° 
‘Equation (39) gives the relationship between the deflections 
for the first and second critical speeds and any assumptions made 
Is must satisfy this relationship. 
Equation (39) may be written 


N 
ZW. 91292 = 0 ° (40) 

; ° 

number of positions considered, 

w.dx = weight of short lengths of the shaft. 


where N 
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The values of W and y, are known from the calculations for 
the first critical speed. Let an assumed deflection curve be drawn 
for the second critical speed in conformity with Fig. 17. 

Let z = distance of the node from the fixed origin. 
y = deflection at any point of the assumed curve. 


The portion of the curve over the distance z is left unchanged 
and the ordinates of the remainder multiplied by a constant K 
chosen that the relationship (40) is satisfied. The constant K is 
found as follows : 


From equation (40), 
BW. e = 2 Wey t+ KZ Wey = 0 
ZW yyy 
ZEW oy» 
where the values of W, y, and y are read for the same points. 


Hence K = (41), 


The corrected curve is then taken as a basis for calculation 
and a value for the second critical speed is computed by formula (20). 


In order to determine the minimum value of the second critical 
speed, it is necessary to make at least three assumptions regarding 
the shape of the elastic curve. The three values of ¢ thus calcu- 
lated are plotted as a function of the distances z of the node points 
from the fixed origin. A smooth curve drawn through the three 
points will yield a minimum value of the second critical speed. 


The procedure may be summarised as follows :— 


(a) Draw the static deflection curve and calculate the funda- 
mental critical speed by formula (20). 


(5) Choose a point on the shaft and assume this is the position 
of the node. In this connection, experience alone is the guide 
to a suitable choice and it may occur that the point chosen will be 
unsuitable and the work will have to be repeated. 


(c) Give the loads on each side of the node point signs in 
conformity with Fig. 17, and draw the deflection curve. In general, 
the position of the node on this derived deflection curve will not 
coincide with that first chosen, but if the two points are in the 
vicinity of each other, the derived deflection curve will be suitable 
for our purpose. 


(d) Calculate the value of the constant K by formula (41), 
and hence correct the ordinates of the assumed: deflection curve. 
It shotild be noted that if the corrected ordinates of the assumed 
curve are such as to produce a very decided discontinuity in the 
corrected curve the work will have to be repeated and a more 
acceptable position of the node chosen. 
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(e) Calculate a value of the second critical speed by formula (20). 

(f) In the same manner, assume two more deflection curves 
and ‘calculate the values of tlie second critical speed for each. It 
may be noted that the choice of the nodes here will be governed 
to a great extent by the results of the first assumption. 

(g) Measure the distance of the nodes on the assumed de- 
flection curves from some fixed origin and plot them against the 
corresponding critical values. The minimum value of the second 
critical speed is then read from the graph, which is drawn smoothly 
through these three points. 
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Production of Accurate Screw Threads, oo 
High Pressure Hydraulic Control Equipment, 
The Fundamentals of Metallurgy, on oo 
Mechanical Methods of Calculation, oo 
Turbo-Supercharging of Internal Combustion Fagines, 
Furnace Brickwork, tts eats wae 
Electrical Protective Equipment, 
Concrete Columns and Founds, wi 
Methods of Approximate Integration, ‘a Pret 
Design Examples for Young Draughtsmen (Reprint), 
Electronic Heating, fen wah sa sag re 
Fundamentals of Gearing (Reprint), ... 
Introduction to Design of Shell Boilers, 
Land Surveying and Levelling, os, 
The Critical or Whirling Speed of Shafts, 


G. H. BrapBury 
H. E, Hutrer 
GEoRGE W. McArp 
B, FELDBAUER 

J. McHarpy Youne 
A. BELL 

G. S. ELFER 


+ C. C. BARNES 
. Frank H. ABRAHAMS 
... RaLpy M. EGGLEston 


WiiiiamM W. CaLpER 
J. McHarpy Younc 
F. G. SPREADBURY 
K. E. HARPER 

W. Scorr WiLson 
M. R. Morton 


GeorGE D. Movat 
T. H. Carr 


F, H. ABRAHAMS 
W. Scott WiLson 


. Duncan CARMICHAEL 


KENNETH J. Fretp 
T. WATERHOUSE 


J. W. K. Murcu 


R. Harrxorr & J. RopGER 
Kennetu G. KErLING 
SypNnEy HarTLey-SMITH 


MackEeNnzIE-KENNEDY 


H. D. Henprick 
T. H. Carr 

W. Scott Wirson 
Gzorce D. Movat 
J. W. Hamitton 
Gro. R. Cooper 
Tuos. C. F. Srorr 
A. GooDALL 

Ian RoBB 
Wi11AM A. NoBLE 


All not otherwise marked are 1/- to members, 2/- others, post free. 
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List of A.E.S.D. Data Sheets. 


Safe Load on Machine-Cut Spur Gears. 

Deflection of Shafts and Beams 

Deflection of Shafts and Beams (Instruction Sheet) 

Steam Radiation Heating Chart. 

Horse-Power of Leather Belts, etc. 

Automobile Brakes (Axle Brakes) 

‘Automobile Brakes (Transmission Brakes) \ Connected. 

Capacities of Bucket Elevators. 

Valley Angle Chart for Hoppers and Chutes. 

Shafts up to 54-in. diameter, subjected to Twisting and Combined 
Bending and Twisting. 

Shafts, 53 to 26-inch diameter, subjected to Twisting and Combined 
Bending and Twisting. - 

Ship Derrick Booms. 

Spiral Springs (Diameter of Rd. or Sq. Wire). 

Spiral Springs (Compression). 

Automobile Clutches (Cone Clutches). 
- i (Plate Clutches). 

Coil Friction for Belts, etc. 

Internal Expanding Brakes, Self-Balancing Brake 
Shoes (Force Diagram) 

Internal Expanding Brakes. Angular Proportions for 
Self-Balancing 

Referred Mean Pressure Cut-Off, etc. 

Particulars for Balata Belt Drives. 

" Square Duralumin Tubes as Struts. 


\ Connected. 


Connected. 


Fi 

1 ” » » ” 

3” Sq. Steel Tubes as Struts (30 ton yield). 

” » (30 ton yield). 

Te ian ” » (30 ton yield). 

vw ” » (40 ton yield). 

Fn ” » (40 ton yield). 

an in » (40 ton yield). 

Moments of Inertia of Built-up Sections (Tables) 2) 


Moments of Inertia of Built-up Sections (Instructions +} Connected. 
and Examples). 
Reinforced Concrete Slabs (Line Chart) 
Reinforced Concrete Slabs (Instructions and Examples) 
Capacity and Speed Chart for Troughed Band Conveyors. 
Screw Propeller Design (Sheet 1, Diameter Chart) 
} Connected 


Connected. 


» (Sheet 2, Pitch Chart) 

” » ” (Sheet 3, Notes and Examples) 

Open Coil Conical Springs. 

Close Coil Conical Springs. 

Trajectory Described by Belt Conveyors. 

Metric Equivalents. 

Useful Conversion Factors. 

Torsion of Non-Circular Shafts. 

Railway Vehicles on Curves. 

Chart of R.S. Angle Purlins. 

Coned Plate Development. 

Solution of Triangles (Sheet 1, Right Angles). 

Solution of Triangles (Sheet 2, Oblique Angles). 

Relation between Length, Linear Movement and Angular Movement 
of Lever. (Diagram and Notes). 
» » ” 5a ae be dt (Chart). 


Helix Angle and Efficiency of Screws and Worms. 
Approximate Radius of Gyration of Various Sections. 
Helical Spring Graphs (Round Wire) 
(Round Wire) Connected. 

a + » (Square Wire) 
Relative Values of Welds to Rivets. 
Ration of Length/Depth of Girders for Stiffness. . 
Graphs for Strength of Rectangular Flat Plates of Uniform Thickness, 
Graphs for Deflection of Rectangular Flat Plates of Uniform Thickness. 
Moment of Resistance of Reinforced Concrete Beams. 
Deflection of Leaf Spring. 
Strength of Leaf Spring. 
Chart showing Relationship of Various Hardness Tests. 
Shaft Horse Power and Proportions of Worm Gears. 
Ring with Uniform Internal Load (Tangential Strain) Z 

it 


Connected 


Ring with Uniform Internal Load (Tangential Stress) 
ress at Bore 


Hub Pressed on to Steel Shaft. (Maximum Tangential 
of Hub). 


68. Hub Pressed on to Steel Shaft. (Radial Gripping Pressure between 


Hub and Shaft). 
Rotating Disc (Steel) Tangential Strain 


= * . rs Stress Connected. 
Ring witb Uniform External Load, Tangential Strain 
» » ” ” ” Stress Connected. 
Viscosity Temperature Chart for Converting Commercial 
to Absolute Viscosities, Connected 


Journal Friction on Bearings, 
Ring Oil Bearings, 


76. Shearing and Bearing Values for High Tensile Structural - 


98-99-100. Boiler Safety Valves. 


Steel Shop Rivets, in accordance with B.S.S. No. 

548/1934, Connected. 
Permissible Compressive Stresses for High Tensile Struc- 

tural Steel, manufactured in accordance with B.S.S. 

548/1934. 
. Velocity of Flow in Pipes for a Given Delivery Cc ted 
79. Delivery of Water in Pipes for a Given Head SANE CHEM 
Working Loads in Mild Steel Pillar Shafts. 
Involute Toothed Gearing Chart. 
. Steam Pipe Design. Chart showing Flow of Steam through Pipes. 
83. Variation of Suction Lift and Temperature for Centrifugal Pumps. 
Nomograph for Uniformily Distributed Loads on British 

Standard Beams, 


” ” ” ” ” 


” ” ” ” ” Connected. 
Notes ee Design and on Use of Data Sheets, Nos. 


Curve Relating Natural Frequency and Deflection, 
Vibration Transmissibility Curve for Elastic Suspension ¢ Connected. 
Instructions and Examples in the Use of Data Sheets, 
Nos. 89 and 90. 
. Pressure on Sides of Bunker. 
93-4-5-6-7. Rolled Steel Sections. 


101. Nomograph Chart for Working Stresses in Mild Steel Columns. 
(Data Sheets are 2d to Members, 4d to others, post free) 
Orders for Pamphelts and Data Sheets to be sent to the Editor, The 


Draughtsman, cheques and orders being crossed ‘“‘Draughtsman 
Publishing Co., Ltd.” 


